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INTRODUCTION 

A set of computer programs has been developed for evaluating the 
performance of powered-lift STOL aircraft. Included are a static per- 
formance summary and dynamic calculations of takeoff and landing per- 
formance. This report describes the input, output, options, and calcu- 
lations for each program. The programs are written in FORTRAN IV and 
are currently available on TSS 360. 

This report is in three independent sections corresponding to the 
three programs: 

(1) Static Performance 

(2) Takeoff Performance 

(3) Landing Performance 

The static performance program computes sets of longitudinal equilibrium 
trimmed flight conditions and displays them in a convenient format. The 
takeoff program estimates the takeoff and climbout maneuver of an aircraft. 
It includes the effects of rotation technique, engine failure, emergency 
thrust, gear retraction and ground effects. The landing program estimates 
the flare and landing roll maneuver. It may also be used to estimate the 
braking distance portion of a rejected takeoff. 

The programs are structured for powered-lift aircraft where the aero- 
dynamics are functions of thrust and velocity as well as angle of attack 
and flap deflection. The engine model in each case will accommodate a 
thrust split between the fan air which interacts with the aerodynamics and 
a vectorable gas generator stream which contributes independently to the 
aircraft accelerations. The engine out capability of these programs is in 
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general restricted to those airplanes where an engine failure and compensa- 
tion produce insignificant lift and drag changes other than those directly 
related to the reduced thrust. The application of these programs to dif- 
ferent powered-lift concepts is discussed in the text. The programs are 
not limited to a specific speed or lift coefficient range. 
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NOTATION 

2 

Horizontal acceleration, body or earth axis, ft/sec 

. . 2 

Vertical acceleration, body or earth axis, ft/sec 

Body station, increasing aft, ft 

Mean aerodynamic chord, ft 

Drag coefficient 

Jet thrust coefficient, Tc/qS 

Lift coefficient 

Lift coefficient at stall 

Pitching moment coefficient 

- / 2 

Acceleration due to gravity, ft/sec 
Altitude 

Moment of inertia in pitch about aircraft C.G., slug-ft 

Incidence of the horizontal tail, deg or rad 

Incidence of the wing, deg or rad 

Pitching moment, ft lb 
Mass of the aircraft, slugs 
Load factor capability, g's 

Number of engines 

Dynamic pressure, psf 

Range, ft 

c 2 

Wing area, ft 

2 

Horizontal tail area, ft 
Specific fuel consumption 

Thrust (cold) interacting with aerodynamics 

Thrust (hot) vectorable 

Flight path velocity, ft /sec, kts 
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Rotation velocity, kts 

Horizontal component of velocity, ft/sec 
Vertical component of velocity, ft/sec 
Velocity at initiation of engine failure, kts 
Gross weight of airplane, lbs 
Water line, increasing up, ft 
Horizontal force, body and stability axes 
Vertical force, body and stability axes 
Angle of attack, fuselage and wing, deg 

Elevator deflection, positive with trailing edge down, deg 

(Augmentor) flap deflection, positive with trailing edge down, deg 

Flight path angle, deg or rad 

Braking coefficient 

Rolling friction coefficient 

Thrust deflection angle, positive down, deg 

3 

Air density, slugs/ft 
Pitch attitude, deg 
Downwash angle, deg 




- 5 - . 

STOL STATIC PERFORMANCE 

The computer program described in this section provides a static 
performance summary for a powered— lift aircraft, where the aerodynamic 
lift, drag, and moment coefficients are functions of thrust coefficient 
as well as angle of attack and flap deflection and where a portion of the 
engine thrust can also be deflected. A longitudinal equilibrium flight 
condition may be specified by angle of attack, flight path angle, thrust 
level and deflection, velocity, and flap deflection. With four of these 
flight variables prescribed by the user, the program varies the other two 
and elevator deflection in an iterative procedure until the equilibrium 
point is reached. The program includes the computation of stall speed 
and maneuver margin and the effects of temperature and altitude. 

The program is an outgrowth of one set up specifically for the Aug- 
mentor Wing Research Aircraft early in 1971. That program was generated 
from the C8-A simulation program by S. M. Sinclair at Ames and then sub- 
stantially revised by The Boeing Company. 1 When the need arose for a more 
generally applicable powered-lift trim program, some further modifications 
were made. In the current program there is a choice between using tail-on 
aerodynamic data or using tail-off data and providing a downwash routine. 
The number of engines can be varied. Although the program was derived for 
an augmentor wing, where the aerodynamics are related to the relatively 
cold thrust from the fan and there is a deflectable hot thrust stream, it 
is applicable to other powered-lift systems, as described in the program 
input section. 


1. The Boeing Company program number TEA-282, documented in their letter 
6-7240-00-152 from R. H. Ashleman to D. D. Few, February 18, 1971. 
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The sections which follow describe the trim algorithm, program input, 
program output, and the programmed equations. 


Trim Algorithm and Stall Speed Calculation 

The trim algorithm seeks a longitudinal equilibrium flight condition 
by iteratively reducing the body x- and z-axis accelerations of the air- 
craft below prescribed tolerance levels. The method is to vary two of the 


flight conditions, alpha, gamma, thrust, and thrust deflection, to reduce 
the accelerations, while varying elevator deflection to maintain zero 
pitching moment. At each stage in the iteration toward trim, the required 
change in acceleration is equal to minus the residual acceleration of the 


previous 

nated 


stage. If the two variables of the trim mode in use are desig- 
and V 2 then the required change may be approximated: 


-a = Aa 
x x 


8a 8a 

8V^ AV 1 + 8vf AV 2 


8a 


8a 


-a = Aa 
z z 


w: 4V i + w 4V 2 


where 


a ,a are the residual accelerations from the previous stage 
x z 

Aa ,Aa are the required acceleration changes to trim 
x z 

V 1 ,V 2 are the trimming variables, two of a, y, v, T^ 

8a 8a are the sensitivities of the accelerations to changes 


x z 

8V . ’ 8V . 

l l 


in and V 2 


These equations are solved for the adjustments to and V 2 to be applied 
before the residual accelerations are recalculated. The actual steps 
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applied are somewhat smaller than these calculated values to avoid con- 
vergence problems due to nonlinearities in the system. 

When the trim mode chosen is to vary flight path angle (or rate of 
climb), the program computes the stall speed in addition to the requested 
equilibrium conditions. The stall speed is found by a Newton-Raphson 
search for the velocity at which the load factor capability, n, is l.Og. 

Expressions for the partial derivatives and the accelerations from 
which they are taken are shown below. (Only partials with respect to the 
two variables of the trim mode in use are required.) The rest of the pro- 
gram equations are included in a separate equations section. 
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Program Input 

The program input has been set up to facilitate parametric variations. 
The major part of the airplane data is input to the program through block 
data routines to be described later. These routines include airplane con- 
stants and all the aerodynamic and engine data in table lookup format. 

Then for any run, only a few data cards are required to select configura- 
tions, velocities, and trim mode, as described below. 


A tri mm ed or equilibrium flight condition is specified by the follow 
ing variables: angle of attack, flight path angle, thrust, thrust deflec- 

tion angle, flight speed and flap deflection angle. During the trimming 
procedure two of these quantities are varied and the remaining four are 
held fixed at their input values. Elevator deflection is varied simultan- 
eously to maintain zero pitching moment. In each of the three trim modes 
available, flap deflection and flight speed are held constant and angle of 

attack is varied. The modes are: 

Mode 1. Thrust and thrust deflection are prescribed. Angle 
of attack and flight path angle are varied to trim. 


Mode 2. Thrust deflection and flight path angle are prescribed. 

Angle of attack and thrust level are varied to trim. 

Mode 3. Thrust and flight path angle are prescribed. Angle of 
attack and thrust deflection are varied to trim. 
Equilibrium conditions are determined for all combinations of the prescribed 
variables. When the trim mode requested is to vary flight path angle, the 
stall velocity corresponding to C Lmax 


is also computed. 



- 9 - 


The input for a run consists of a title card and a NAMELIST data set 
(INPUTS) which must conform to the 360 NAMELIST convention. The user need 


only set 

those variables for which he wishes to change 

from 

the default 

values . 

The most frequently used input variables are: 



"INPUTS" 

MAX. 

DEFAULT 

NAME 

DESCRIPTION 

NO. 

VALUE 

DFLAP 

Flap settings, deg. 

10 

0 . 

PWR 

Power settings, in whatever units used 

10 

0 . 


in block data II 



GAMZ 

Flight path angles, deg. 

10 

0 . 

R0C 

Rates of Climb, ft /min 

10 

0 . 

DEFL 

Hot thrust (Th) deflection angles, deg. 

10 

0 . 

VKN0TS 

Flight path velocities, knts. 

10 

0 . 

HZ 

Altitude, ft. 

1 

0 . 

WT 

Gross weight, lbs 

1 

As set in 




block data IV 

DELTAT 

Std day temp deviation, °F 

1 

0 . 

C0NR0C 

Constant r/c code 

1 

0 . 


= 0. GAMZ input 




= 1. ROC input 



M0DTRM 

Trim condition code (Integer) 

1 

-2 


< 0 Vary a & y 

(Set 

PWR, DEFL) 


= 0 Vary a & hot thrust 

(Set 

GAMZ or ROC, DEFL) 


> 0 Vary a & deflection angle 

(Set PWR and GAMZ or ROC) 

I0E0 

Engine out code (Integer) 

1 

0 


= 0 All engines operating 




= 1 One engine out (for propulsion systems 

where insignificant 


asymmetry is produced, see page 12) 

Additional variables may be changed in namelist INPUTS from their values 

I *' 

set in the block data routines. These variables, as defined in the block 
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data sections, are: 

CDT, CDE, CLDE, IW, IT, AOT , AIT, SW, SHT, CBAR 

Data sets may be stacked, but note that input values are carried 
over between input data sets. 

Once the block data routines have been compiled, a run deck will be 
of the following form: 

LOAD block data routines 
CALL program 
Title card 
^INPUTS 

Set desired variables name = value separated by commas 
$END 

Additional title and input sets 
%END 

Aerodynamic Data 

The program is most simply used if the aerodynamic data is available 
with tail on. If the data is not available in that form, then a downwash 
calculation must be provided by the user. Only a non-functioning dummy 
routine exists in the program, because the subroutine originally in the 
program was specific to the Augmentor Wing Research Aircraft. The user 
must provide a subroutine of the following form: 

SUBROUTINE DWNWSH 
C OMMON / AER0XX / EP S , XA , Z A 

COMMON / TRMC0N / D F , NUR , CNU , SNU , THROTL , GAMMA , ALPHA , 

1C A, S A, THETA, ALT , DELT AT , DE , MODTRM , GAIN, G , OMASS 
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Any statements required to compute EPS = epsilon, the 
downwash angle, at the horizontal tail, in degrees. 

RETURN 

END 

Other labelled commons which may be needed appear in the block data 
description. 

Effect of the Propulsion System 

The program was derived for an augment or-wing configuration where the 
aerodynamics are related to the relatively cold thrust from the engine fan. 
This thrust is labelled T . The aerodynamic coefficients are functions of 
jet thrust coefficient C , as well as angle of attack and flap deflection, 

J 

where C = T /qS. The thrust from the hot core gases, which is labelled 
J c 

T h> can be deflected through the angle v, where zero v is the undeflected 

case. The hot thrust T^ makes a separate contribution to the normal and 

longitudinal accelerations of the aircraft. 

The program can be used for other powered-lift aircraft, or even for 

conventional airplanes, if one treats these other propulsion systems as 

special cases of the augmentor-wing situation. In externally blown flap 

configurations, the aerodynamics are functions of total thrust. To use 

this program set T to the total mixed thrust and T, to zero. Then the 

c n. 

calculated thrust coefficient will be (Total Thrust) /qS and there will be 
no deflected thrust. In this case, the user must select the power setting; 
and alpha, gamma, and elevator angle will be computed for equilibrium. 

For configurations where the aerodynamics are unrelated to the pro- 
pulsion system, use T^ as the total thrust and set T c to zero. In this 
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/ 


case the total thrust can be deflected through the angle v. The thrust 
coefficient is zero. All three trim modes are available in this situ- 
ation. A conventional airplane would correspond to this second situation 
with no thrust deflection. The following sketch illustrates the different 
applications that can be made. These are applicable to the takeoff and 
landing performance programs as well. 

The engine out capability of these programs is in general restricted 
to those airplanes where an engine failure and compensation produce insig- 
nificant lift and drag changes other than those directly related to the 
reduced thrust. In cases where a significant asymmetry is produced, a 
separate set of aerodynamic data must be provided for the engine out sit- 
uation, including the required corrections for lateral and directional 


trim. 
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Output 

A set of trim conditions corresponding to the requested velocities is 
generated for each combination of the flaps, powers, flight path angles, 
and/or thrust deflections input. A summary of the output symbols is shown 
below. 


Symbol Dimension 

Run Title Information: 


Wt 



GW 


LBS 


ALT 


FT* 


DELTA T 


°F 


FLAPS 


DEG 


I TAIL 


DEG 


A0T 


RAD 


I WING 




S WING 


FT* 


SH TAIL 


FT 


CBAR 


FT 


FROM C.G. 

AFT 

DOWN 


MOM CTR 

XX 

XX 

CBAR 

TAIL 

XX 

XX 

CBAR 

ENG THR 

XX 

XX 

CBAR 

ENG INL 

XX 

XX 

CBAR 


Output Head Line: 


V 

a 

Y 


v 



£ 


pwr 



c 


VEL Knots TAS 

GAMMA 

DEG 

R/C 

Ft /Min 

NU 

DEG 

ALPHA F 

DEG 

THETA 

DEG 

DE 

DEG 

ALPHA T YdLwash?!5 

EPS i used 


CJ 

2 

AX 

FT/SEC„ 

AZ 

FT/ SEC 

POWER SETTING 


HOT THRUST 

LBS 

COLD THRUST 

LBS 


Description 


Gross weight 
Altitude 

At in reference to a standard day 
Flap setting 

Tail incidence in ref. to W.L. 

Zero lift line of the tail 
Wing incidence in ref. to W.L. 

Reference wing area 
Horizontal tail area 
Mean aerodynamic chord 

Positions with respect to the center of 
gravity of moment center, horizontal tail, 
engine thrust point, and engine inlet point 


Flight path velocity 
Flight path angle 
Rate of climb 

Hot thrust deflection angle 
Fuselage angle of attack 
Body attitude 
Elevator deflection angle 

Horizontal tail angle of attack 
Downwash angle at MAC horizontal tail 
Cold thrust coefficient 
Residual tangential acceleration 
Residual normal acceleration 
Power setting in units used in input 
Hot thrust 
Cold thrust 
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Symbol 


Dimension 


Description 



RAM DRAG 
CL TRIM 


LBS 

L wb 

CLWB N 

i 

1 


D wb 

CDWB 

j 

r 


Kb 

CMWB J 




CL M TRIM 
CL M WB 




ALPHA M 


DEG 

Lmax 




N 


G 


V-VMIN 
V/V MIN 


knots TAS 


ENGINE LIMIT 



Ram drag 

Trimmed lift coefficient 

With downwash calculation: 

tail-off lift, drag, and moment coeff. 

Without downwash calculation: 
tail-on untrimmed coefficients 

Trimmed max. lift coefficient 

Tail off (or tail-on untrimmed) max. lift coeff. 
Fuselage angle of attack at C Lmax 

Load factor capability 

Difference between "VEL" and stall speed 

Ratio of "VEL" and stall speed 

(Blank) - No limit encountered 

MAX - Power setting limit exceeded 

TBLE - Power setting outside table range 
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BLOCK DATA ROUTINES 


A description of the block data subroutines required to input aero- 
dynamic and engine data and airplane constants is given on the following 
pages. Each routine is of the form: 

BLOCK DATA 
COMMON statements 

DATA statements to set all variables described 
END 


Block Data I. Aerodynamic Data 
Common Statements : 

COMMON /AR0TB1/NAR0 (3) , ATAB (31) , CLTAB (900) , CDTAB (900) ,CMTAB (900) 

C OMMON / ATBL S / CLT AB 1(150), CLTAB 2 (150) ,CLTAB3(150) ,CLTAB4(150) , 
1CLTAB5(150) ,CDTAB1(150) ,CDTAB2(150) ,CDTAB3(150) ,CDTAB4(150) , 

2CDTAB5 (150) .CMTAB1 (150) , CMTAB2 (150) , CMTAB3 (150) , CMTAB4 (150) , 
3CMTAB5(150) ,CLTAB6(150) ,CDTAB6(150) ,CMTAB6(150) 

C0MM0N/CLMX/NCLX(2) ,TABCLX(21) ,CLWBMX(90) ,NAX(2) ,TABAX(21) ,ALPHMX(90) 
COMMON / TLON / T AILON 


Data to be Input : 

Stability axis aerodynamic coefficients 

C L* C D ,C M = ftn <V C J’ 6 P C J = V qS ^ See section on Propulsion system) 

a W = °F + \ 


^Lmax* “C 
TAILON 


Lmax 


= ftn(Cj, < 5 f ) 

= 1. to input tail-on data 

= 0. tail- off (In this case see section on aerodynamic data.) 
Number of values of at which C^> C^, to be input (>2) 
Number of values of Cj. at which C L> Cp, C M to be input (>2) 


NAR0(1) 

NAR0(2) 
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NAR0(3) Number of values of 6 f at which C L> C D ,C M to be input 

(2<NARO(3) <6) 

NAR0(1) + NAR0(2) + NARO(3)<31 
ATAB 


NAR0(1) * NAR0(2)<150 


CLTABl 

CDTAB1 

C11TAB1 


NAR0(1) values of u , NAR0C2) values of C Jt and MAR0C3) . 

values of 6 f , each set monotonically increasing (<31 entries; 

Matrix of lift coefficients C L = ftn^.Cj) for the first 
6^, varying most rapidly 

Matrix of drag coefficients C D = ftnC^.Cj) for the first 

S,., a rT varying most rapidly 
f W 

Matrix of pitching moment coefficients- C M = ftn(a w> C J ) for 
the first S f , varying most rapidly 

CLTAB2-6 CDTAB2-6, CMTAB2-6 Aerodynamic coefficients at the remaining 

flap settings, parallel to the above. Use as many as required. 

These data are combined into the vectors in common / AROTBL/ and three- 
dimensional lookup carried out (TVIN) . 

NCLX(l) Number of values of Cj at which C Lmax to be input (>2) 

NCLX(2) 

NCLX(l) + NCLX(2) <21 NCLX(l) * NCLX(2)<90 
TABCLX 


Number of values of 6 at which C Lmax to be input (>2) 


CLWBMX 


NCLX(l) values of Cj followed by NCLX(2) values of S f , each 
set monotonically increasing 

Matrix of values of C Ljmx - ftnfCj.dj). C_, varying moat rapidly 


MAX, TABAX Same as NCLX, TABCLX above for at - ^ Lmax 

Lmax 


ALPHMX Matrix of values of a 

rapidly 


- ftn(C ,6_) , C T varying most 
C T j r j 

Lmax 


Block Data II Thrust 
Common Statements ; 

COMMON /PROT AB/NPR0(3) ,PTAB(22) ,HOTTH(300) ,COLDTH(300) , 
1EMFTAB(300) ,ETATAB(300) ,BL0SS(12) ,NOE,PMAX 
COIiMON / BOUND /VI , VG , V2 
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NOE 

NPRO(l) 

NPRO(2) 

NPRO(3) 


Data to be Input : 

T c , T h Ramdrag - ftn(power setting, altitude, velocity) 

(Be sure you have read section on effect of propulsion system.) 

T = T * eta * bloss * engines operating 
c C i 

T = T * bloss * engines operating 
h 

Ramdrag = Ramdrag * engines operating 

Number of engines (at full power) 

Number of power settings at which thrusts to be input 
(2<NPR0(1)<12) 

Number of altitudes at which thrusts to be input (>2) 

Number of velocities at which thrusts to be input (>2) 

NPRO(l) + NPR0(2) + NPRO(3)i22 NPR0(1) * BPR0C2) * m*0O)<J00 

S(3) values °of 
set monotonically increasing 

Maximum allowable value of power setting 

„ c t — ffnfnnwer alt, vel) for one engine, in lbs, 
power X varying _ most 'rapidly ? 1 fhen anitud-fhen velocity 
(TVIN 3-dimensional lookup) 

Matrix of T c = ftn (power, alt, vel) parallel to T h 

Matrix of Ramdrag >= ftn(power,alt,vel) parallel to T h 

Matrix of cold thrust efficiency = ftn (power, alt, vel) parallel 
to T h 

Bleed loss efficiency - ftn( P ower setting), NPR0(1) entries 
Applied to both hot and cold thrust 


PMAX 

HOTTH 

COLDTH 

EMFTAB 

ETATAB 

BLOSS 


Bounds for stall velocity search. 
VI 


a f* -/ope for which aerodynamic data is avail 

S3 r^locity 

search ihould be bounded.) 



V2 


Corresponding maximum speed 
A starting velocity within the above range 


VG 


Block Data III . Moment Arms 
Common Statement : 

COMMON/ CGTRMS/CGT 1, CGT2, CGT4, CGT5 ,CGT6,CGT 7, XE,ZE 

Data to be Input : 

Distances from C.G. to various engine positions, tail and moment center 
in terms of mean aerodynamic chord. 


WL - water line, increases up 


BS = Body station, increases aft 


CGT1 

CGT2 

CGT4 

CGT5 

CGT6 

CGT7 

XE 

ZE 


(WLCG 

(BSCG 

(BSCG 

(WLCG 

(WLCG 

(BSCG 

(BSCG 

(WLCG 


WLMC) / C 
BSMC) /c 

BSTQ) / c 
WLTQ) /c 

VJLEI) / c 
BSEI) / c 

BSENG) /c 
WLENG) /c 


Moment center position 

f 

Horizontal tail quarter chord position 
Engine inlet position 

Engine hot thrust exit position or point 
of application of hot thrust vector 


Note that on the program output these distances are shown as how far 


down, "how far" aft, so that for the horizontal positions the signs will be 


shown opposite to those input above. 


Block Data IV . Airplane Constants and Initializations 
Common Statements : 

COMMON/ APCON / IT , 1$ , SW , SHT , STOSW, CBAR, WT , AOT , AlT , ALPHT 
REAL IT,IW 

COMMON/COEFF/ CLWB , CDWB , CMWB , PDCL ( 3) , PDCD (3) , T CLCOF , TCDCOF , CDT , CDE , 


1CLT.CLDE 
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COMMON/NOZANL/BETNOZ .PHINOZ ,PSINOZ 

COMMON /TRMCON /DF , NUR , CNU , SNU , THROTL , GAMMA , ALPHA , C A , S A , THETA , ALT , 

DELT AT , DE , MODTRM , GAIN , G , OMAS S 

COMMON/TRMPRM/DFLAPS(10) ,PWR(10) ,GAMZ(10) , ROC (10) ,DEFL(10) 

1VKNOTS ( 10 ) , CONROC , HZ 

COMMON /ENGN/TH, TC , CJ ,XENG , ZENG ,MENG , IOEO , RDRAG , C5X, C5Z 

Required Initializations : 

DATA HZ, MODTRM, DELTAT.G/O. ,-2,0. ,32.174/ 

Default altitude to zero, trim mode to varying gamma, 
standard-day temperature deviation to zero, and set g 

DATA DFLAP , PWR , DEFL , GAMZ , ROC , VKNOTS / 10 *0 . , 10*0 . , 10*0 . , 10*0 . , 10*0 . , 10*0 . / 

Default the input vectors to zero 

DATA CONROC/O. / 

Default constant rate of climb flag to zero 


DATA IOEO/O / 


Default to no engine out 


Data to be Input : 


Airplane dimensions and derivatives: 

WT Default value for gross weight of airplane 

. 2 

SW Wing area, ft 


SHT 

CBAR 

IW 

CDT 

CDE 


Horizontal tail area 

Mean aerodynamic chord, ft 

Incidence of the wing, radians 

Contribution of tail to total drag coefficient in terms of 
wing area. (Zero for tail— on aero data.) 

Contribution of elevator to total drag coefficient in terms 
of wing area 
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CLDE 

IT 

BETNOZ 

PHINOZ 

PSINOZ 

Following 

AOT 


3C / 36 , per radian, in terms of tail area 
L e 

Incidence of the horizontal tail, radians. (Will be printed 
but only used with downwash calculations.) 

Installation offset angles of the hot thrust deflection nozzle 
radians. (See equations.) 


quantities only used with downwash calculation. 

Zero lift angle of attack of the horizontal tail, radians 

Linearized lift curve slope of the horizontal tail, per radian 
in terms of tail area 


AIT 



- 26 - 


STOL TAKEOFF PERFORMANCE 

TAKOFF is a computer program for estimating the takeoff and climbout 
maneuver of a STOL aircraft with powered-lift characteristics. The pro- 
gram includes the effects of rotation technique, engine failure, emergency 
thrust, gear retraction, and ground effects. It is an outgrowth of a sim- 
ilar unpublished program for conventional aircraft by V. R. Corsiglia of 
Ames Research Center 

For the powered-lift aircraft, aerodynamic lift, drag, and moment are 
functions of thrust level and velocity as well as angle of attack and flap 
deflection. Tail-on aerodynamic coefficients C^, C^, are input to the 
program in tabular form as functions of thrust coefficient and angle of 
attack for a specified flap deflection. The engine model has. been set up 
to accommodate the augmentor wing system, where the total thrust is split 
between a hot thrust part which is exhausted through vectorable nozzles 
and a cold thrust part which interacts with the aircraft aerodynamics. 

The program is applicable to other powered— lift aircraft, however, as 
described in the input section. 

For purposes of the computer analysis, the takeoff maneuver is divided 

into four segments: 

(1) ground roll to rotation velocity 

(2) rotation to the commanded maximum angle of attack, usually 
including liftoff 

(3) transition to the climbout condition 

(4) climbout 

The options available and calculations carried out will be described sep- 
arately for each segment. The equations for all segments are included in 
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the computations section. 

Ground Roll to Rotation Speed 

The first segment of the takeoff maneuver, to be called ground roll 
for convenience, extends from brake release until the aircraft has accel- 
erated to the selected rotation velocity, V R . The distance is integrated 
from an initial time, position, and velocity, which need not be brake re- 
lease. This allows the calculations to be restarted, as demonstrated in 
the sample output. When the velocity, V^, at which an engine fails is 
reached, thrust adjustment begins. Thrust and ramdrag from the failed 
engine are reduced to zero over the winddown interval, DTI. Meanwhile, 
after a delay, DT2, for the pilot to respond to the failure, the remain- 
ing engines are moved to full emergency thrust over a time period, DT3. 

This period of thrust change, shown schematically in figure 1, normally 
extends into the rotation segment. 

Prior to rotation, angle of attack is held at a constant input value 
and elevator deflection is taken to be zero. Ground effect factors as a 
function of altitude are applied to and from the initiation of ground 
roll until the aircraft is out of ground effect. 

Rotation 

When the aircraft reaches the rotation velocity, V R , an elevator step 
change 6e 0 is input. Provided this elevator step is sufficient to trim 
the aircraft pitching moment and still lift the nose gear at V R , TAKOFF 
iterates to find the elevator pattern for which the aircraft will just 
rotate to the commanded maximum angle of attack. The quantity determined 



Normal 

power 



Other Engines 



t yl t yl + DT2 tpre + DT3 

(tpre) 


'hen an engine fails . 
g) 




Technique for Determining 


slope = yj 



a. Initial guess t^ = t^ 

b. Compute rotation 


+ DTR1 


c. If at t - ^ > .03 rad/sec 

decrease t^. 

(if t 1 < t RQT reduce 6e 0 ) 
da „ , 

If < 0 before a - a . 

L dt max- rot’ 

i 

increase t. a 


d. Repeat b and c until 

\^\ < .03 rad/sec at t = t^ 


Figure 2. Time history of & & and trial and error technique 
for determining t, . 


fiZ 
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by the program is the time period the full elevator step is to be held 
before the elevator is returned at a specified rate. The procedure is 
shown in figure 2. Liftoff usually occurs near the end of rotation. 

Transition 

During the transition segment, the aircraft moves from the maximum 
rotation position to the climbout condition. Ground effects can be in- 
cluded until the end of transition. Angle of attack is commanded by the 
user as a function of altitude during this segment. When the gear retrac- 
tion height is reached, TAKOFF begins to retract the gear by removing the 
gear increment from the drag coefficient. After a time period, DT4, the 
gear is fully retracted. During transition and climbout, elevator deflec- 
tion is that which is required to maintain zero pitching moment. 

Climbout 

Several modes of climbout are available: constant rate of climb, 

constant pitch attitude, constant load factor, constant flight path angle, 
and constant velocity. In each case TAKOFF adjusts angle of attack as re- 
quired to maintain the climbout mode as closely as possible without losing 
velocity. Optionally the power setting may be reduced for climbout. There 
will be discontinuities between transition and climbout. 

Program Input 

Input to the program is in two sections. The aerodynamic coefficients 
and engine data are input through a block data subroutine described in the 
next section. Airplane constants and the variables to specify a particular 
takeoff situation are input on a series of twelve data cards as shown on 


the following pages . 
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Card 

1 

2 


Input Data Set 


Column 

Format 

Variable 

Use 

1-80 

20A4 

IDENT 

Title information 

1-8 

F8. 

SMAX 

Terminal distance, ft. 

9-16 

F8. 

W 

Airplane gross weight, lbs. 

17-24 

F8. 

SREF 

2 

Reference area for C L »C D ,C M in ft 

25-32 

F8. 

VR 

Speed at which rotation begins, knots 

33-40 

F8. 

VI 

Speed at which engine failure begins, 
knots (Enter high value if engine 
failure is to be omitted) 

41-48 

F8. 

ALMAX 

a , maximum angle of attack reached 

.max’ . . 
in rotation 

49-56 

F8. 

SFC 

Specific fuel consumption (lb of 
fuel/hr)/(lb of total thrust) 

57-64 

F8. 

HMAX 

Terminal altitude above runway, ft 

65-72 

F8. 

PWR 

Power setting during normal operation 

73-80 

F8. 

EPWR 

Power setting during emergency oper- 
ation after engine failure 

1-8 

F8. 

NU 

Angle of thrust vector from airplane 
reference line, deg. Positive down 

9-16 

F8. 

ALO 

Angle of attack during ground roll 
prior to time VR reached, deg 

17-24 

F8. 

HGEAR 

Altitude above runway at which to 
begin gear retraction, ft. Set = 0 
for fixed gear 

25-32 

F8 . 

DCDG 

AC D , Cd increment due to gear 

bei§f a Sown, a positive quantity. 
Set = 0. for fixed gear 

33-40 

F8. 

DT4 

Time required for gear to retract, sec 

41-48 

F8. 

CBAR 

Mean aerodynamic chord, ft 
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Card 

Column 

Format 

Variable 

Use 

3 

49-56 

F8. 

ICG 

Moment of inertia i| pitch about air- 
plane C.G., slug-ft 


57-64 

F8. 

XCG 

Dimensionless distance with respect to 
CBAR, from moment center aft to C.G. 


65-72 

F8. 

XGR 

Dimensionless distance w.r.t. CBAR 
from moment center aft to main gear 


73-80 

F8. 

ZGR 

Dimensionless distance w.r.t. CBAR 
from moment center down to main gear 

4 

1-8 

F8. 

ZCG 

Dimensionless distance w.r.t. CBAR 
from moment center down to C.G. 


9-16 

F8. 

XEI 

Dimensionless distance w.r.t. CBAR 
from C.G. aft to engine inlet 


17-24 

F8. 

ZEI 

Dimensionless distance w.r.t. CBAR 
from C.G. down to engine inlet 

5 

1-8 

F8. 

DTI 

Time from initial engine failure 
until complete loss of thrust, sec 
(>0) 


9-16 

F8. 

DT2 

Time from initial engine failure 
until pilot has reacted and completed 
throttle adjustments for emergency 
thrust, sec (>0) 


17-24 

F8. 

DT3 

Time required by engines to reach 
full emergency thrust after throttles 
adjusted, sec (>0) 


25-32 

F8 . 

MU 

Rolling friction coefficient of gear 
on the runway 


33-40 

F8. 

CMDE 

Derivative of C^ 3C^/35_, 1/deg, in 
terms of dimensionless distance 


41-48 

F8. 

CLDE 

Derivative of C^ 3C^/36 e , 1/deg 


49-56 

F8. 

CDDE 

Derivative of C^ 30^/36^, 1/^e 8 


57-64 

F8. 

DEO 

Value of 6. after step change at 
beginning of rotation, deg. 


65-72 

F8. 

DET 

36 /3t, deg/sec 
e 


73-80 

F8. 

NOE 

Number of engines 
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Card 

Column 

Format 

Variable 

Use 

6 

1-5 

15 

NOH 

Number of altitude entries in tran- 
sition tables for ALGE, FCLGE, FCDGE 
( 2 <NOH< 7 ) 


6-10 

15 

MXPRNT 

= 1 to print trial rotation calcula- 
tions. Otherwise not 


11-15 

15 

TCPRNT 

= 1 to print cold thrust, T £ . De- 
faults to print hot thrust, 

7 

1-8 

F8. 

XTH 

Dimensionless distance w.r.t. CBAR 
from C.G. aft to the engine hot 
thrust reaction point 


9-16 

F8. 

ZTH 

Corresponding distance down from C.G. 


17-24 

F8. 

DTRl 

Initial guess of length of time 
(Tl - T ) that total elevator step 

will be neld, seconds. Defaults to 
1. second 


25-32 

F8. 

IW 

Incidence of the wing. 


33-40 

F8. 

VZ 

Initial velocity, kts (VZ<VR) 


41-48 

F8. 

HZ 

Runway elevation, ft 


49-56 

F8. 

AUGRAT 

Augmentation ratio. (Used in extra- 
polating C , C beyond table limits - 
see computations section. Set AUGRAT=0. 
for no extrapolation.) 


57-64 

F8. 

SZ 

Initial distance, ft 


65-72 

F8. 

TZ 

Initial time, sec 

Transition: 




8 

1-70 

7F10 . 

H 

NOH values of altitude above runway 
at which transition tables input, in 
ascending order and no two identical. 
The aircraft is out of ground effect 
and climbout begins at H(NOH). 
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Card 

Column 

Format 

Variable 

Use 

9 

1-70 

7F10 . 

ALGE 

NOH values of angle of attack as 
functions of altitude during transi- 





tion 

10 

1-70 

7F10 . 

FCLGE 

NOH corresponding multiplicative cor- 
rections to C L for ground effect 

11 

1-70 

7F10 . 

FCDGE 

NOH corresponding multiplicative cor- 
rections to C D for ground effect 


The factors FCLGE and FCDGE are 
during ground roll and rotation. 

also applied to C^ and C^ 

Climbout: 




12 

1-6 

A6 

MODE 

Designated mode of climbout. Values 
ROC, THETA, LDFCTR, GAMMA, VEL cor- 
respond to constant rate of climb, 
pitch attitude, load factor, flight 
path angle, and velocity, respectively 


7-12 

F6. 

VALUE 

The value of the mode of climbout 
with units ft/min (ROC), degrees 
(THETA, GAMMA) , knots (VEL) 


13-16 

FA. 2 

Thrust 

Factor 

Multiplicative factor to reduce 
power setting for climbout. (Thrust 
is discontinuous from transition to 


climbout . ) 

Sets of data cards may be stacked. Follow the last set with: %END 


Effect of the Propulsion System 

The program was set up for an augmentor-wing configuration where the 
aerodynamics are related to the relatively cold thrust from the engine fan 
This thrust is labelled T q . The aerodynamic coefficients are functions of 
jet thrust coefficient C Jt as well as angle of attack and flap deflection, 
where C = T /qS. The thrust from the hot core gases, which is labelled 
T , can be deflected through the angle v, where zero v is the undeflected 
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case. The hot thrust T h makes a separate contribution to the normal and 

longitudinal accelerations of the aircraft. 

The program can be used for other powered-lift aircraft, or even for 
conventional airplanes, if one treats these other propnlsion systems as 
special cases of the augmentor-wing situation. In externally blow, flap 
configurations, the aerodynamics are functions of total thrust. To use 
this program set T^ to the total mixed thrust and T h to zero. Then the 
calculated thrust coefficient will be (Total Thrust) /qS and there will 
be no deflected thrust. 

For configurations where the aerodynamics are unrelated to the pro- 
pulsion system, use T h as the total thrust and set T c to zero. In this 
case the total thrust can be deflected through the angle v. The thrust 
coefficient C T is zero. A conventional airplane would correspond to this 

second situation with no thrust deflection. 


Limitations 

The engine out capability of the takeoff program and the others is in 
general restricted to those airplanes where an engine failure and compen- 
sation produce insignificant lift and drag changes other than those directly 
related to the reduced thrust. In cases where a significant asymmetry is 
produced, a separate set of aerodynamic data must be provided for the en- 
gine-out situation, including the required corrections for lateral and 
directional trim. For the takeoff calculation, all-engine and engine-out 
results would have to be combined to estimate takeoff performance with 
failure at a particular speed. This would only be approximate. 
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Further limitations on the class of airplanes for which the programs 
are appropriate are minimum speed and rotation pattern. For the takeoff 
and landing calculations, extrapolation of aerodynamic data to high Cj 
values is provided to allow computation during ground roll. Nominal zero 
velocity is taken at one f.p.s. in calculating these coefficients. The 
takeoff calculation is restricted to that class of airplanes for which 
rotation for takeoff can be represented as a control input 6, and esti- 
mates of andC Lj are available. While this does not require 6, 

to he elevator deflection, the control input will be referred to as ele- 
vator (6 ) throughout the report. 

Engine and Aero dynamic Data 

A block data subroutine is required to input the engine thrust and 
aerodynamic tables. It must include the following statements: 

BLOCK DATA 

COMMON /AERO/ TH1TAB(150), THCTAB(150), XTAB(25), 

1XYZ(3) , NXYZ (3) , PARDER(3) , NALCJ(2), ALCJ(25), XY(2), 

2CLTAB(150), CDTAB(150), CLOSS(150), BLOSS(IO), NPWR, 

3N0TH, NOV, NAL, NCJ, RAMDR(150), DFLAP, CMTAB(150) 

c C D , and C M ar<? defined as functions of alpha^ and Cj, where Cj = 

(total cold thrust) / (Q * S ref ) . tail-on untrimmed. To define C L and 

C for a given flap setting, the following variables must be set: 

D 
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NAL Number of values of alpha^^ (where = a p + i y ) for aerodynamic 

tables C L , C D , C M = ftn (a^Cj) . Must be >2. 

NCJ Number of values of Cj for aerodynamic tables. Must be >2. 

(NAL + NCJ < 25 and NAL * NCJ < 150) 

ALCJ NAL values of alpha (in deg) followed by NCJ values of C^, each 
set being in increasing order with no two identical. 

CLTAB Table of C L values as function ( a w » C j)> varying most rapidly. 

CDTAB Table of C^ values as function (a^,Cj), varying most rapidly. 

Gear down. 

CMTAB Table of CM values as function ( a w > C j) parallel to C^.C^. 

DFLAP Flap setting for which C^, C^ input, deg. Used in estimation of 

C , C when C out of range. See computations section. 

Li D J 

The thrust provided by each engine is defined as a function of power, eleva- 
tion, and velocity, modified by loss factors. For each engine. 

T, = THl(pwr, H, V) * BLOSS(pwr) 
h 

T = THC(pwr, H, V) * CL0SS(pwr, H, V) * BLOSS(pwr) 
c 

Ramdrag = RAMDR(pwr, H, V) 

The following variables must be set: 

NPWR Number of power settings for thrust tables. Must be >2 and <10. 

NOTH N umb er of elevation settings for thrust tables. Must be >2. 

NOV Number of velocity settings for thrust tables. Must be >2. 

(NPWR + NOTH + NOV <25, NPWR * NOTH * NOV < 150) 

XTAB NPWR values of power followed by NOTH values of elevation above 

sea level (ft), followed by NOV values of velocity (knots). Each 
set must be in increasing order with no two identical. 
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TH1TAB Table of hot thrust values for one engine as function (pwr,H,V), 
power varying most rapidly, then elevation, then velocity. 

THCTAB Parallel table of cold thrust values for one engine. Total cold 
thrust is used to compute Cj. 

CLOSS Parallel table of multiplicative correction factors for losses 
in cold thrust. 

BLOSS NPWR multiplicative correction factors for bleed valve loss as a 
function of power setting. Applied to both hot and cold thrust. 

RAMDR Table of ramdrag levels for each engine parallel to the thrust 
tables . 


Program Output 

Output from’ TAKOFF consists of a summary of the input data and 
history. The variables which are computed and printed at intervals 
second are: 


velocity 

dV 

dt 


thrust 


gamma 


alphay 
C D 


distance 
da 


dt 

C, 


elevation 

d 2 a 

dt 2 

weight 


a time 
of . 1 


6 load factor 

e 

In addition the time history shows the beginning of engine failure and 

rotation, liftoff, and the beginning of climbout. The thrust level printed 

is the hot thrust % unless the flag TCPRNT has been set to request T c to 

n 

be printed. 

Optionally the output includes the repeated time histories calculated 
during the iteration of the rotation sequence. This is rarely required. 
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Computations 

The following integrations are carried out by the program: 
V = 7^7 flight speed 

y = flight path angle 



H - M 
H “ J dt 


range 

elevation 


a 


F 


da 

dt 


da _ ,d_a 
dt J dt 


angle of attack 

time rate of change of angle of attack 


The required derivatives at any time T are computed by the equations 
shown below. The expressions common to all segments of the takeoff maneu- 
ver are shown before those specific to certain segments. 

Fuel consumption 


W = W . , - Sfc * (T, +T 
old h c 

Dynamic pressure 

p = ftn(H) 

q = .5pV 2 


T - T 
3600 


old 


Thrust 

T^ = ftn(pwr,H,V) * bleed loss factor * N 

T c = f tn(pwr,H,V) * bleed loss factor * cold loss factor * N 
Ramdrag = ftn(pwr,H,V) * N 

(Thrusts are adjusted at engine failure as shown in figure 1.) 
Cj • T c /qS 
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Trim elevator 

6 _ = [C qS (xcg cosct - zcg sim p ) + C D qS(xcg sinotp+ zcg cosctp) 


trim 


F 

+ C M qS + T^(zth cosv - xth sin v) - Ram(zei cosotp+ xei sinoip)]/ 
9C. 


< - jr ^ 

e 


On ground: 


3C. 


M 


6 = 6 + Ntfrc(xgr - xcg + p(zgr - zcg))/( - qS) 

trim trim e 


where Ntfrc = C^qS + sin(Op + v) — 
Aerodynamic adjustments 
Gear drag reduction 


W 


T - T 

FDCDG = - ( £ 


) AC n 

gear time Dgear 


for T> T and limited to - AC 

g Dgear 


Ground effect factors 


FCL, FCD = ftn(H) 

For climbout FCL, FCD = 1. 

Extrapolation for C T beyond table limit C 

J J I 

AC l = Aug. Ratio * (Cj - Cj )sin (a y + <5 f ) 
AC d = -Aug. Ratio * (Cj - Cj )cos(a w + <S f ) 


Aerodynamic coefficients 

8C 

C = [ftn(a C )+ AC l + ST 6 e ] * FCL 

extrap e 

aC D 

C_ = [ftn(a C )+ AC d + FDCDG + — 6 £ ] * FCD 

w extrap e 

C M = ftn( V C J } 

Net force 

Ntfrc = C L qS + T h sinCa^ + v) - W 
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Derivatives 


ip = & [T h cos (ot p + v) - C D qS - Wsiny - Ramdrag] 
On ground: 

QL = 4 ^ + y * Ntfrc * g/W 
dt dt 


^ = g[Ntfrc/W + 1 . - cosy]/V 


On ground: 

£-<>• 


ds „ 

— = Vcosy 
dt 


dH 

dt 


= Vsiny 


Quantities dependent on takeoff segment 


Ground roll (V < V R ) 


6 e - 0. 


a F = a. 


da 

dt 


d 2 a 

dt 2 


= 0 . 


Rotation (a p < a rotation ) 


6 e = 6 e step input 


3(5 


If T > T x 6 e = 6 e<> + - 9 — (T - T x ) 

where T, is end of held elevator step 

3C m - 

Moment = e “ & e ) 35 ^ Sc 
trim e 

,2 

— j = Moment/ I 

dt 2 g 2 

Nose gear lifts if 77 = is positive. 

at dt > 
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Transition (H < H , 

max ground effect 


6=6 

e e 


trim 
otp = ftn(H) 


da 

dt 


d 2 a 

dt 2 


= 0 . 


Climbout (H > H , 

- max ground effect 

Pwr = Pwr * thrust factor for climbout 

Determine alpha according to climbout mode, updating 

C T , C at each iteration 

Lj U 

3C, 


C L = ftn < a W’ C J )+ 36 6 e . 

e trim 

SC D 

C D ' £tn( V C J >+ FDCDG + 35“ 


da 

dt 


d 2 a 

dt 2 


= 0 


'trim 


dV 


Reduce a if necessary to maintain > 0. 

Climbout modes 

Constant velocity 

dV 

Search for a such that = 0. 

Constant ROC 

Eq = (CpC.osy - C D siny)qS/W + T h s.in(a F +v+y) /W - .1 - Ramdrag/W * 

Search for a such that Eq = 0. 

Constant pitch attitude 

a - 0 - y where 0 is the specified value 
F 

Constant load factor 

Load factor = (C L qS+ T h sin(a F +v) ) / (Wcosy) 

Search for a such that Load factor = specified value 


siny 
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Constant gamma 



dy 

Search for a such that — = 0. 



Summary 

of Moment Arm Definitions 




Moment center aft to C.G. 

Bs 

b cg 

Bs 

s mc 

XCG 

c 



Moment center down to C.G. 

Wl 

L mc 

Wl 

eg 

ZCG 

c 


XGR 

Moment center aft to main gear 

Be - 

b gJ 

c 

Bs 

°mc 


Moment center down to main gear 

i 

% 

3 

W T 

V 

ZGR 

c 

pt 

Bs - Be 

b eng b cg 

C.G. aft to engine hot thrust reaction 

XTH 

c 

ZTH 

C.G. down to engine hot thrust reaction pt 

Wt — Wt 

‘-■eg ‘-'eng 

c 


C.G. aft to engine -inlet point 

Bs el - 

BS 

eg 

XEI 

c 



C.G. down to engine inlet point 

wl 

eg 

WL ei 

ZEI 

c 



where Wl = water line, increases up 

Bs = body station, increases aft 
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STOL LANDING PERFORMANCE 

LAND is a simple performance program with trimming, landing and open- 
loop longitudinal plane maneuvering capabilities. The program is specif- 
ically structured for a class of power-augmented lift aircraft for which 
the aerodynamic lift and drag coefficients are functions of thrust level 
and velocity as well as angle of attack and flap deflection. The mathe- 
matical model has only two longitudinal degrees of freedom with the pitch- 
ing moment equation being replaced by a command angle of attack during the 
maneuvering flight computations. Tail-on lift, drag, and moment data 
must be tabulated in coefficient form as functions of angle of attack, 
thrust coefficient, and flap deflection. The thrust coefficient expresses 
the influence of engine thrust level on the aerodynamic effectiveness of 
the lifting surfaces. The engine model will accommodate a split in the 
total thrust between a hot thrust part which is exhausted through vector- 
able nozzles and a cold thrust part which interacts with the aircraft 
aerodynamics. It is applicable to other powered- lift systems as is de- 
scribed in the input section. 

The program performs the iterative computation to trim the aircraft 
in an equilibrium flight condition prescribed by program input. At touch- 
down, the program executes a landing roll computation which may include 
braking, thrust reversal, thrust change and lift dumping, all prescribed 
by program input. Maneuvering flight segments are controlled by subroutine 
CNTRL. The supplied CNTRL maneuver is a simple flare control. It allows 
for increments in angle of attack and thrust specified by the user as a 
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function of time. For any other flight maneuver the user may supply a 
subroutine specifying angle of attack as a function of time, height, range 
or any other flight variable(s) accessible to the computation. 

The program description which follows includes an outline of the 
program input requirements, a discussion of the trimming technique, and 
description of the landing roll calculation and of the control subroutine. 


Trim Algorithm 

A trimmed or equilibrium flight condition is specified by: angle of 

attack, flight path angle, thrust, thrust deflection, flight speed, and 
flap deflection angle. During the trimming procedure, two of these var- 
iables and elevator are varied until the residual x— and z— axis accelera- 
tions of the aircraft have been reduced below prescribed tolerance levels. 
In the three programmed trim modes, angle of attack is varied along with 
gamma, thrust, or thrust deflection. 

At each stage in the iteration toward trim, the required change in 
acceleration is equal to minus the residual acceleration of the previous 
stage. If the two variables of the trim mode in use are designated 
and V 2 then the required change may be approximated: 


-a 


x 


= Aa 


x 


3a 

3V^ 1V 1 + 


9 a 

3V 2 AV 2 


-a 


= Aa 


9 a 

9^ AV 1 + 


9a 

9V7 AV 2 


where 


a . a are the residual accelerations from the previous stage 
x z 

Aa . Aa are the required acceleration changes to trim 
x z 
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V 

8a 
x 

8V . * 
1 


are the trimming variables, two of ct, y» v > 

8a are the sensitivities of the accelerations to changes 
z_ 

8V. . „ , „ 

l in and V ^ 


These equations are solved for the adjustments to and to be applied 
before the residual accelerations are recalculated. The actual steps 
applied are somewhat smaller than these calculated values to avoid con- 


vergence problems due to nonlinearities in the system. 

Expressions for the accelerations in stability axes and the required 
partial derivatives are shown in the equations list. (Only partials with 
respect to the two variables of the trim mode in use are required.) 


Landing Calculation 

Xn the normal mode of operation the aircraft is first trimmed in the 
equilibrium flight condition specified and flown through a maneuver pre- 
scribed by the control (CNTRL) subroutine. If the resulting flight path 
reaches ground level, a landing roll and braking sequence is initiated. 
The prescribed rolling friction coefficient takes effect immediately upon 
touchdown. The braking coefficient may be applied after a specified de- 
lay interval. Change in power setting, thrust reversal and lift dumping, 
if desired, begin after a reaction or delay time from touchdown and are 
completed within a time interval prescribed by input. The user specifies 
the final power setting and hot thrust vector angle to be used during 
braking. The flag STODUMP controls the thrust coefficient after touch- 
down, maintaining the appropriate value for the power setting (NODUMP-1) 
or reducing the coefficient to zero, simulating spoiling of the power 


induced lift (N0DUMP=0). 
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The trimming and landing roll functions may be used independently by 
setting flag NOLAND (see input list) . 

The equations used are shown in the equations section. 

LAND Input 

Input to the program is in two sections. The engine and aerodynamic 
characteristics of the airplane are input in tabular form through a block 
data subroutine. Library routines BVIN and TVIN are used for two- and 
three-dimensional table lookup. Characteristics of a particular landing 
are input by a set of five control cards. Control sets may be stacked. 


Input Data Set 


Card 

Column 

Format 

Variable 

Use 

1 

1-80 

20A4 

TITLE 

Any desired title 

2 

1-8 

F8. 

ZOBS 

Obstacle height 


9-16 

F8. 

MUR 

y = rolling friction coefficient 
R 


17-24 

F8. 

MUB 

y = braking coefficient 
B 


25-32 

F8. 

DTB 

Delay time after touchdown before 
brake applied 


33-40 

F8. 

W 

Aircraft weight 


41-48 

F8. 

S 

Reference area for aerodynamics 


49-56 

F8. 

DT 

Time step for integration 


57-61 

15 

NOLAND 

= 1 Trim only 
= 0 Regular trim and land 
= -1 Omit trim and calculate brak- 
ing distance (ZOBS should be I 

3 

1-8 

F8. 

VAPP 

Approach speed, knots 
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Card Column 

Format 

Variable 

Use 

9-16 

F8. 

APPPWR 

Approach power setting. If MODTRM = 0 
(see card 4) , this value of APPPWR 
will be replaced by trim value 

Thrust reversal 

control 



17-24 

F8. 

REVPWR 

Power setting at completion of thrust 
reversal for braking 

25-32 

F8. 

DTREV 

Delay time after touchdown to begin- 
ning of thrust reversal, sec 

33-40 

F8. 

REVTIM 

Time interval for change in power, 
lift dumping, and deflection of thrust 
to reverse position, sec 

41-48 

F8. 

REVNU 

Deflection of hot thrust at full 
reversal, deg 

49-53 

15 

NODUMP 

= 1 to prevent lift dumping 
= 0 to dump lift by sloping cold 

thrust (and thereby C^) to zero 
at thrust reversal 

54-58 

15 

NOREV 

= 1 hot thrust not reversed. Nu 
maintained at nu trim. (Dump 
and power change continued) 

= 0 to reverse thrust (deflection 
to REVNU) 

Trim control 




4 1-5 

15 

MODTRM 

= -1 to trim by varying alpha and 
gamma 

= 0 to trim by varying alpha and 
hot thrust, 

= 1 to trim by varying alpha and 
thrust deflection v. 

6-13 

F8. 

GAMMA 

Fixed value of gamma for trim of 
MODTRM = 0 or 1 

14-21 

F8. 

NU 

Fixed value of thrust deflection if 
MODTRM = -1 or 0 

22-29 

F8. 

DFLAP 

Flap deflection 

30-34 

15 

NOE 

Number of engines 



- 49 - 


Card 

Column 

Format 

Variable 

Use 

Flare 

control 

using provided subroutine CNTRL 

5 

1-8 

F8. 

FLRTIM 

Length of time over which Aa 
for flare is applied, secs 


9-16 

F8. 

ALPINC 

Aa for flare, deg 


17-24 

F8. 

DTFLAR 

Delay from T^before begin flare, secs 


25-32 

F8. 

PWRTIM 

Length of time over which power incre- 
ment for flare applied, sec 


33-40 

F8. 

PWRINC 

Power increment for flare, units of 
power setting 


41-48 

F8. 

DTPWR 

Delay from T before flare power 
change begun, sec 


Follow the last stacked data set with: %END 


Effect of the Propulsion System 

The program was derived for an augmen tor-wing configuration where the 
aerodynamics are related to the relatively cold thrust from the engine fan. 
This thrust is labelled T^ The aerodynamic coefficients are functions 
of jet thrust coefficient C Jf as well as angle of attack and flap deflec- 
tion, where C = T /qS. The thrust from the hot core gases, which is 
labelled T h> can be deflected through the angle v, where zero v is the 
undeflected case. The hot thrust \ makes a separate contribution to the 
normal and longitudinal accelerations of the aircraft. 

The program can be used for other powered-lift aircraft if one treats 
these other propulsion systems as special cases of the augmentor-wing sit- 
uation. In externally blown flap configurations, the aerodynamics are 
functions of total thrust. To use this program set T c to the total mixed 

thrust and T, to zero. Then the calculated thrust coefficient will be 
h 
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(Total Thrust) /qS and there will be no deflected thrust. In this case, 
the user must select the power setting, and equilibrium alpha, gamma, 
and elevator angle will be computed. (MODTRM = -1) 

For configurations where the aerodynamics are unrelated to the pro- 
pulsion system, use T^ as the total thrust and set T^ to zero. In this 
case the total thrust can be deflected through the angle v. The thrust 
coefficient C^. is zero. All three trim modes are available in this sit- 
uation. 


Airplane Description 


Engine and aerodynamic data are input in tabular form through a block 
data subroutine. Some geometric constants are also included in that 
routine. The subroutine must be of the following form: 


BLOCK DATA 

COMMON/C/CLTAB(400),CDTAB(400) ,NALCJ (3) ,ALCJ (25) , 

ID FLAP , AUGRAT 

COMMON /THR/THTAB (100) ,TCTAB(100) ,RAMTAB (100) ,NPV(2) , 

1PV(25) ,NOE,CLOSS (100) ,BLOSS(10) 

COMMON /M/ CMTAB (400 ) , DETRIM , XCG , ZCG , XTH , ZTH , 

1CMDE , CLDE , CDDE , XEI , ZEI , IW 
REAL IW 

Data statements to set all input variables to be described. 


END 

The input variables to be included are defined below. Be sure you have 

read the section on the effect of the propulsion system before continuing. 

Aerodynamics (Tail on data, untrimmed, C^ = T c /qS) 

NALCJ(l) Number of values of alphay at which C L > C D » C M 

to be input. (ct w = oi^, + i w ) (>2) 
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NALCJ(2) Number of values of C at which C ,C ,C to 

be input. (>2) J 

NALCJ (3) Number of values of 6 f at which C L ,C D ,C M to 

be input. (>i2) 

NALCJ (1) * NALCJ (2) * NALCJ (3) < 400, 

NALCJ (1) + NALCJ (2) + NALCJ (3) < 25 


ALCJ 


CLTAB 

CDTAB 

CMTAB 

AUGRAT 


CLDE 

CDDE 

CMDE 

Geometry 

XCG 

ZCG 

XTH 

ZTH 


NALCJ (1) values of a (deg), followed by 
NALCJ(2) values of C , followed by NALCJ(3) 
values of 6^ (deg) at which C^,Cp,C^ to be 
input, each set monotonically increasing 

Matrix of C^ as a function of a w , C , 6 , 
varying most rapidly, then C^, then 

Parallel matrix of C^ values, gear down 

Parallel matrix of C_, values 

M 

Augmentation ratio, used to extrapolate C^, 

C D for C beyond table range. The method of 
extrapolation is shown in the equations section. 
If no augmentation, input 1. 

To avoid extrapolation of C ^>C D input 0. 

9C /36 , 1/deg 

Li 6 

3C D /96 e , 1/deg 

3C^/96 e , in terms of wing area, 1/deg 


Distance from C.G. 
to moment center of 
data 


Distance from C.G. 
to engine hot thrust 
exit point 




XEI 
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ZEI 


Engine 

NOE 

NPV(l) 

NPV(2) 

PV 

THTAB 

TCTAB 

RAMTAB 

CLOSS 

BLOSS 

T, = NOE 
h 

T = NOE 
c 

Note that 



Distance from C.G. 
to engine inlet point 


- Waterline, positive up 
Bg - Body station, positive aft 

Number of engines 

Number of power settings at which thrust to 
be input. (>2) 


Number of velocity settings at which thrust 
to be input. (>2) 


NPV(l) values of power level followed by 
NPV(2) values of velocity (knots) at which 
thrust to be input, each set monotonically 
increasing 

Matrix of hot thrust in pounds for one engine 
as a function of power and velocity, power 
varying most rapidly 


Parallel matrix of cold thrust values 


Parallel matrix of ramdrag values 


Parallel matrix of multiplicative cold loss 
factors 

Vector of multiplicative bleed loss factors 

as function of power setting, applied to 

both T, and T 
h c 

THTAB (pwr,vel) * BLOSS (pwr) 

' TCTAB (pwr, vel) * CLOSS (pwr, vel) * BLOSS(pwr) 

NPV(l) * NPV(2) <100, NPV(l) + NPV (2) <25, NPV(1)<10, 
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The Control Subroutine 

The control subroutine provided allows for flare control by increments 
in angle of attack and thrust as a function of time. The timing and incre- 
ments are input by the user (see data card 5). This routine may be replaced 
by a user subroutine of the following form: 


SUBROUTINE CNTRL 
COMMON /CNT/ ALTRIM, NUTRIM, 
1PWRINC , DTPWR 

COMMON /D/ TDWN, TH, TC, MU, G, 


FLRTIM, ALPINC, DTFLAR, PWRTIM, 
Q, OMASS, LIFT, DRAG, VEC(51), 


1RAMDR, DTREV, NODUMP 

COMMON /T/ VA, GM, ALPHA, NUR, APPPWR, PWR, MODE, S, W, DF, NOCONV 


REAL NUR, NUTRIM 

Any desired manipulation of 
based on altitude (VEC(6)), 
touchdown time (TDWN) etc. 


ALPHA (radians) or PWR setting 
velocity (VA fps) , time (VEC(2)), 


RETURN 

END 

If no control statements are inserted, the trim value of alpha will 
be maintained. (ALTRIM (radians)) For a braking distance calculation 
without trim, zero alpha will be maintained. 
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APPENDIX 

Sample output from the three programs is shown on the following 
pages. The output corresponds to the baseline configurations of a jet- 
powered STOL augmentor wing design for which an analysis of takeoff and 
landing performance was previously reported. All calculations for that 
analysis were done with the computer programs described in this report. 
The sample output in this appendix shows: 

(1) Static performance - engine out full power in the take- 
off configuration, 41.1°6^ and thrust undeflected. 

(2) Takeoff performance - nominal all engine takeoff to 35 
foot barrier (initial portion of ground roll omitted). 

(3) Landing performance - nominal landing from 35 foot barrier 
at 70. 6°6 f and T/W = .2, hot thrust deflected 90°. 

The aerodynamic data and engine assumptions used for these runs are 
available in the referenced report. 


1. Post, S. E., Gambucci, B. J., and Holzhauser, C. A., An Analysis of the 
Takeoff and Landing Performance of a Jet-Powered STOL Augmentor Wing 
Design , NASA TM X-62,176, August 1972. 



£ INPUTS 

DFL AP = 41.09999, 9*0.0 
PWR* 38.0, 9*0.0 
GAMZ= 10*0.0 
R0C= 10*0.0 
DEFL* 10*0.0 

VKNOTS* 50.0, 60.0, 70.0, 80.0, 90.0, 100.0, 110.0, 3*0.0 

WT = 48000.0 

HZ* 0.0 

DEL TAT* 0.0 

COT* 0.0 

COE* 0.0 

CLDE* 1.695999 

IW= 0.0 

IT* 0.0 

AOT* 0.0 

AIT* 0.0 

SW= 600.0 

SHT* 156.0 

CBAR* 9.679999 

ZSTAR* 0.0 

0R3* 0.0 

0R3CC* 0.0 

IOEO* 1 

CONROC* 0.0 

MOOTRM* -2 

£ENO 



I 

Gn 

00 

I 





RUN TITLE 

- ENGINE OUT 

TRIM AT 

BASE THRUST 

FLAP 

61.1 

GW 

68000. 

LBS 



S 

WING 60 

ALT 

0. 

FT 



FROM C.G. 

DELTA T 

0. 

DEG F 




MOM CTR 

FLAPS 

61.1 

DEG 




TAIL 

I TAIL 

0.000 

DEC- 




ENG THR 

AOT 

0.000 

RAD 




ENG INL 

I WING 

0.000 

DEG 



CLWB IS TA 

VEL 

GAMMA 

R/C 

NU 

ALPHA F 

THETA 

DE 

( KNOTS) 

(DEG) 

( FT/M IN ) 

(DEG) 

(DEG) 

(DEG) 

(DEG) 

53.9 

-5.21 

-695 .6 

Q.O***** 5 **** 

27.56 

-25.70 

60.0 

-0.20 

-21.1 

0.0 

22.13 

21.96 

-16.97 

70.0 

3.57 

661.5 

0.0 

12.59 

16.16 

-5.51 

80.0 

5.61 

792.0 

0.0 

5.82 

11.63 

0.68 

90.0 

6.31 

1001.5 

0.0 

1.70 

8.01 

6.66 

100.0 

6.22 

1097.9 

0.0 

-1.05 

5.17 

7.87 

110.0 

5.81 

1128.5 

0.0 

-3.10 

2.72 

10.37 


600.00 F T **2 


AFT 

- 0.000 

6.080 

0.208 


down 

0.000 C BAR 
-1.180 CBAR 
0.176 CBAR 
0.176 CBAR 


SH TAIL 
CBAR 


156.00 FT**2 
9.68 FT 


TAIL-ON UNTRIMMED 

ALPHA T EPS 
( DFG ) (DEO) 

*«**t**4sc****ti* 
<: *: z*;* — £ $$$ ***** ££ 
£*£££»*£**£*£*££ 
££*££££ **££#£»*£ 
**************** 


CJ AX A l 

FT/SEC2 FT/SEC2 
1. 66 0 ** **************** 


1 .325 
0.978 
0.753 
0.598 
0.686 
0.606 


0.0069 

0.0236 

0.0239 

0.0179 

0.0126 

0.0163 


0.0619 

0.0670 

0.0363 

0.0087 

-0.0233 

-0.0665 


VEL 

POWER 

HOT THRUST 

COLD THRUST 

(KNOTS) 

SETT ING 

(LBS) 

(LBS) 

53.9 

38.0 

3180.0 

9690.0 

60.0 

38.0 

3180.0 

9690.0 

70.0 

38.0 

3195.0 

9738.0 

80.0 

38.0 

3210.0 

9786.0 

90.0 

38.0 

3225.0 

9836.0 

100.0 

38.0 

3260.0 

9882.0 

110.0 

38.0 

3255.0 

9930.0 


RAM DRAG 

CL TRIM 

CD TRIM 

CL WR 

CDWR 

CMWB 

(LBS) 

0 .9 8R7********* 

0.9887 

-0.9563 

1323.00 

6.3910 

0. 2605 

6.5055 

0.2*05 

-0.5866 

1699.60 

6.7359 

-0.1626 

6.7779 

-0.1626 

-0.2358 

1675.80 

3.6650 

-0. 2672 

3.6611 

-0.2672 

-0.0155 

1852.20 

2.89 26 

-0.238R 

2.8567 

-0.2388 

0.1295 

2028.60 

2.3532 

-0.1975 . 

2.2927 

-0.1975 

0.2366 

2205.00 

1.9522 

-0.1563 

1.8726 

-0.1563 

0.3162 


I 

,171 
f VO 

I 


VEL 

(KNOTS) 

CL M TRIM 

CL M WB 

ALPHA M 
(DEG) 

N 

(G) 

V-VMIN 

(KNOTS) 

V/ V M I N 

53.9 

7.820 

8.018 

32.76 

1.002 

0.0 

1 .0000 

60.0 

7.086 

7.273 

32.07 

1.115 

6. 1 

1.1128 

70.0 

6.369 

6.551 

31.03 

1.358 

16. 1 

1 .2982 

80.0 

5.831 

6.003 

30.26 

1.621 

26. 1 

1 .6837 

90.0 

5.391 

5.561 

29.55 

1.893 

36.1 

1.6692 

100.0 

5.031 

5.156 

28.96 

2.175 

66. 1 

i .8566 

110.0 

6.762 

6.871 

28.68 

2.686 

56. 1 

2.0601 



STOL TAKF-OFF SI MULATI OM 
ROT AT I ON WITH STICK POSITION ASSIGNFP 
TRANSITION WITH ASSIGNED CL. CD, ALPHA 
CL IMBDUT WITH ASSIGNED RATE DF CL I MR, 
THFTA, LOAD FACTOR, FLIGHT PATH 
ANGLE, OR VELOCITY 


TAKEOFF 


BASE CASE TAKEOFF T/W=3a% ALL ENGINES 


INPUT UNITS.. .VEL-KNOTS, TIME-SEC, DIST-FT, ANGLE-DEG, PENS 1 TY- SLUGS /FT**3 , T HRUST , WE IGHT-LRS, 


RATE OF CL I MB-ET /MIN 


SMAX 

3000. 

W 

ALMAX 

12.0 

SEC 

ALO 

0.0 

HGEAR 

XCG 

0.000 

XGR 

MU 

0.030 

CMDE 

NOE 

A. 

XE! 

XTH 

0.207 

ZTH 

SEC 

0.000 

AUGRAT 


48000. 

SREF 

600. 

0.000 

HM AX 

40.0 

5. 

DCDG 

0.020 

0.620 

ZGR 

1.240 

-0.02930 

CLOE 

0.0078 

-0.775 

ZEI 

0.176 

0.176 

DTR 1 

1.000 

1.3 

ZCG 

0.000 


VR 

72.0 

VI 

Nil 

0.00 

PWR 

DT 4 

5.00 

CBAR 

DTI 

1.0 

DT2 

code 

0 .00000 

DEO 

IN 

0.0 

VZ 

SZ 

856.7 

TZ 


999.0 


38.0 

E PWR 

38.0 

9.7 

ICG 

0.2300 

1.0 

DT3 

1.0 

1 

o 

DET 

O 

O 

70.0 

13.6 

HZ 

o 

O 


H 

ALGF 

FCLGE 

FCOGE 


0. 

12.0 

1.000 

1.000 


35. 

6.0 

1.000 

1.000 


MODE= VEL 


VALUE= 79.70 THRUST FACTOR 1.00 


CD 

O 

I 



T H X L v. 


TIME 

VEL 

GAMMA 

SEC 

KNOTS 

OEG 

13.6 

70.0 

0.00 

13.7 

70.4 

0.00 

13.fi 

70.9 

0.00 

13. 9 

71.3 

0.00 

14.0 

71 .7 

0.00 

begin rotation 


14.1 

72.0 

0.00 

Tl = 

14.5175 


14.1 

72.1 

0.00 

14.2 

72 . 6 

0.00 

14.3 

73 .0 

0.00 

14.4 

73.4 

0.00 

14.B 

73.8 

0.00 

14.6 

74.2 

0.00 

14.7 

74.6 

0.00 

14.fi 

75.0 

0.00 

14.0 

75.3 

0.00 

15.0 

75.7 

0.00 

15.1 

76.1 

0.00 

15.2 

76.4 

0.00 

LIFTOFF 



15.3 

76.6 

0.00 

15.3 

76.7 

0.00 

15.4 

77.0 

0.06 

15.5 

77.3 

0.19 

15.6 

77.6 

0.39 

15.7 

77.fi 

0.66 

15.fi 

78.0 

0.98 

15.9 

78.2 

1.35 

16.0 

7fi.3 

1.77 

16.1 

7fi .4 

2.14 

16.1 

78.5 

2.19 

16.2 

78.6 

2.54 

16.3 

78.7 

2.89 

16. A 

78.8 

3.23 

16.5 

78.9 

3.58 

1 6.6 

79.0 

3.9 1 

16.7 

79.0 

4.25 

16.fi 

79.1 

4.58 

16.9 

79.2 

4.90 

17.0 

79.2 

5.21 

17.1 

79.3 

5.52 

17.2 

79.3 

5.82 


OIST 

HT 

OV/DT 

FT 

FT 

FT / SEC ** 2 

856.7 

0.0 

7.29 

868.6 

0.0 

7.27 

flfl0.5 

0.0 

7.24 

892.5 

0.0 

7.21 

904.6 

0.0 

7.18 

912.8 

0.0 

7.09 

TR = 14 . 

0675 


916.8 

0.0 

7.08 

929.0 

0.0 

7.03 

941.3 

0.0 

6.98 

953.7 

0.0 

6.91 

9 66.1 

0.0 

6.81 

978.6 

0.0 

6.70 

991.2 

0.0 

6.57 

1003.8 

0.0 

6.41 

1016.5 

0.0 

6.25 

1029.3 

0.0 

6.10 

1042.1 

0.0 

5.92 

1055.0 

0.0 

5.72 

1063.7 

0.0 

5.54 

1067 .9 

0.0 

5.43 

10R0.9 

0.0 

5.07 

1094.0 

0.0 

4.61 

1107.0 

0 . 1 

4.17 

1120.2 

0.2 

3.71 

1133.3 

0.4 

3.19 

1146.5 

0.7 

2.70 

1159.8 

1.0 

2.28 

1170.8 

1.4 

2.01 

1173.0 

1.5 

2.08 

1186.3 

2.0 

1.91 

1199.5 

2.7 

1.74 

1212.fi 

3.4 

1.59 

1226 . 1 

4.2 

1.44 

1239.4 

5.1 

1 . 30 

1252. R 

6.0 

1.17 

1266 . 1 

7.0 

1.06 

1279.4 

8.1 

0.95 

1292.7 

9.3 

0.86 

1306. 1 

10.6 

0.77 

1319.4 

11.9 

0.69 


ALPHA 

OAL/OT 

0 * 92 AL 

DFG 

DFG/SFC 

OEG / SFC**2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0 . 0 

0.0 

0.0 

0 .0 

0.0 

0.0 

0.0 

0.0 

0.0 

5.4 

0.0 

0.2 

5.7 

0.1 

0.8 

6.7 

0.2 

1.5 

7.8 

0.4 

2.4 

9 .1 

0.6 

3.4 

10.7 

1.0 

4.4 

10 .1 

1.5 

5.4 

9.2 

2.1 

6.3 

8.5 

2.8 

7.1 

7.9 

3.5 

7.9 

7.3 

4.4 

ft. 6 

6.9 

5.2 

9.2 

6.5 

5.9 

9.7 

6.2 

6.2 

9.8 

4.8 

7.2 

10.1 

0.7 

8.2 

10.0 

- 3 . 5 

9.2 

9.4 

-7 .8 

10.1 

8.4 

- 12.1 

10.9 

7.0 

-16 .4 

11.5 

5.1 

- 20 . 7 

11.9 

2.9 

- 25.0 

12.0 

0.6 

- 28 . 4 

11.7 

0.0 

0 .0 

11.6 

0.0 

0.0 

11.5 

0.0 

0.0 

11.4 

0.0 

0 . 0 

11.3 

0.0 

0.0 

11. 1 

0.0 

0.0 

11.0 

0.0 

0.0 

10.8 

0.0 

0.0 

10.6 

0.0 

0 .0 

10.4 

0.0 

0.0 

10.2 

0.0 

0 .0 

10.0 

0.0 

0 . 0 


THRUST 

CO 

CL 

LBS 

4260 . 

- 0.905 

3.427 

4261 . 

- 0.891 

3.410 

4262 . 

- 0.876 

3.394 

42 63 . 

- 0.861 

3.378 

4263 . 

- 0 . R47 

3.362 

4264 . 

- 0. 838 

3.001 

4264 . 

- 0.833 

2.996 

4265 . 

- O.filfi 

2 .989 

4266 . 

- 0. 801 

2.992 

4267 . 

- 0.781 

3.008 

4268 . 

- 0.759 

3.039 

4268 . 

- 0.733 

3.117 

4269 . 

- 0.705 

3.215 

4270 . 

- 0.674 

3.329 

4271 . 

- 0.642 

3.453 

4271 . 

- 0.612 

3.584 

4272 . 

- 0 . 580 

3.727 

4273 . 

- 0.544 

3.876 

4273 . 

- 0.516 

3.975 

4273 . 

- 0.502 

4.025 

4274 . 

- 0.458 

4.179 

4275 . 

- 0 . 40P 

4.358 

4275 . 

- 0.368 

4.509 

4276 . 

- 0.329 

4.651 

4276 . 

- 0.288 

4.772 

4276 . 

- 0.254 

4.873 

4277 . 

- 0.232 

4.951 

4277 . 

- 0.223 

4.992 

4277 . 

- 0.235 

4.712 

4277 . 

- 0.237 

4.697 

4277 . 

- 0.240 

4.681 

4278 . 

- 0.244 

4.663 

4278 . 

- 0.249 

4.645 

4278 . 

- 0.254 

4.625 

4278 . 

- 0.261 

4.603 

4278 . 

- 0.268 

4.581 

4278 . 

- 0.277 

4.557 

4278 . 

- 0.286 

4.532 

42 79 . 

- 0.295 

4.506 

4279 . 

- 0.305 

4.478 


W 

OE 

LOAD 

LBS 

OEG 

FACTOR 

48000. 

0.0 

0.713 

48000. 

0.0 

0.718 

48000. 

0.0 

0.723 

48000. 

0.0 

0.729 

48000. 

0.0 

0.734 

48000. 

-45.0 

0.660 

48000. 

-45.0 

0.662 

48000. 

-45.0 

0.668 

48000. 

-45.0 

0.677 

48000. 

-45.0 

0 .688 

48000. 

-45.0 

0.703 

48000. 

-41.3 

0.730 

48000. 

-36.8 

0.761 

48000. 

-32.3 

0.797 

48000. 

-27.8 

0.836 

48000. 

-23.3 

0.877 

48000. 

-18.8 

0.922 

48000. 

-14.3 

0.969 

48000. 

-11.3 

1 .000 

48000. 

-9.R 

1.016 

48000. 

-5.3 

1 .064 

48000. 

-O.P 

1.119 

48000. 

3.7 

1.167 

48000. 

8.2 

1 i 2 1 1 

48000. 

12.7 

1.250 

48000. 

17.2 

1.283 

48000. 

21.7 

1.309 

48000. 

25.5 

1.324 

48000. 

-6.4 

1 .251 

48000. 

-6.3 

1.251 

48000. 

-6.2 

1 .250 

48000. 

-6.1 

1.248 

48000. 

-6.0 

1.267 

48000. 

-5.9 

1.244 

48000. 

-5.8 

1.241 

48000. 

-5.6 

1.237 

48000. 

-5.5 

1.233 

48000. 

-5.6 

1.228 

48000. 

-5.2 

1.223 

48000. 

-5.1 

1.217 


' I 

O' 

i-* 

i 



I A KC- Off F 


17.3 

79.3 

6.12 

1332.7 

13.3 

17.4 

79.4 

6.40 

1346.1 

14.7 

17.5 

79.4 

6.67 

1359.4 

16.3 

17.6 

79.4 

6.93 

1372.7 

17.9 

17.7 

79.5 

7.17 

1386.1 

19.5 

17.8 

79.5 

7.41 

1399.4 

21.2 

17.0 

79.5 

7.63 

1412.7 

23.0 

18.0 

79.5 

7.83 

1426.0 

24.8 

18.1 

79.5 

8.03 

1439.3 

26.6 

18.2 

79.5 

8.20 

1452.6 

28.5 

18.3 

79.6 

8.36 

1465.9 

30.5 

18.4 

79.6 

8.51 

1479.2 

32.4 

18.5 

79.6 

8.63 

1492.5 

34.4 

VELOCITY = 79.7 KNOTS 

35 FT ABOVE RUNWAY 

BEGIN 

CONSTANT 

18.5 

79.6 

8.67 

1496.1 

35.0 

18.6 

79.6 

8.78 

1505.8 

36.5 

18.7 

79.6 

8.95 

1519.1 

38.6 

18.8 

79. 6 

9. 10 

1532.4 

40.7 

40. 

FT ALTITUDE REACHED 


18.8 

79. 6 

9.05 

1528.2 

o 

• 

o 

4- 


0.61 

9.7 

0.0 

0.0 

42 79 

0.53 

9.5 

0.0 

0. 0 

4279, 

0.46 

9.2 

0.0 

0.0 

4279 

0.40 

8.9 

0.0 

0. 0 

4279, 

0.36 

8.7 

0.0 

0.0 

4279 

0.32 

8.4 

0.0 

0.0 

4279, 

0.30 

8.1 

0.0 

0.0 

4279 

0.28 

7.8 

0.0 

0. 0 

4279, 

0.30 

7.4 

0.0 

0.0 

42 79 

0.31 

7.1 

0.0 

0. 0 

4279 

0.33 

6.8 

o.n 

0.0 

42 79 

0.35 

6.4 

0.0 

0.0 

4279 

0.39 

6.1 

0.0 

0 .0 

4279 

VELOCITY 

11 

sO 

• 




0.41 

6.0 

0.0 

0. 0 

4279 

-0.00 

7.1 

0.0 

0 .0 

4279 

0.00 

6.8 

0.0 

0. 0 

4279 

0.00 

6.6 

0.0 

0.0 

42 79 

O 

o 

• 

0 

1 

6.6 

© 

O 

c 

c 

42 79 


-0.314 

4.446 

48000. 

-4.9 

1.210 

-0.322 

4.414 

48000. 

-4.8 

1.202 

-0.332 

4.380 

48000. 

-4.6 

1.194 

-0.341 

4.345 

48000. 

-4.5 

1.186 

-0.351 

4.308 

48000. 

-4.3 

1.177 

-0.362 

4.271 

48000. 

-4.2 

1.167 

-0.373 

4.233 

48000. 

-4.0 

1.157 

-0.384 

4.193 

48000. 

-3.8 

1.147 

-0.398 

4.143 

48000. 

-3.7 

1.134 

-0.410 

4.101 

48000. 

-3.5 

1.123 

-0.422 

4.059 

48000. 

-3.4 

1.112 

-0.434 

4.016 

48000. 

-3.2 

1.101 

-0.447 

3.972 

48000. 

-3.0 

1 .089 


-0.450 

3.960 

48000. 

-2.9 

1.086 

-0.411 

4.097 

48000. 

-3.5 

1.125 

-0.422 

4.062 

48000. 

-3.4 

1.116 

-0.432 

4.030 

48000. 

-3.2 

1.107 

-0.428 

4.042 

48000. 

-3.3 

1.110 


W= 48000. LBS 


P-H 0 = 


0.002375 SLUGS/FT**3 



landing 

T/W=.20 

70.6 FLAP 

68 KT 





N) $ 



WT = 48 000. 


MUR= 0 

.030 

APPPWR= 

20.0 

VAPP 

= 

68 . 

OFLAP = 70.6 




S = 600. 


MU B= 0 

.300 

REVPWR= 

20.0 

RE VT 1 M= 

l.on 

RFVNIJ = 114.0 




DT = 0.2 


DTB = 

1.0 

DTREV= 

0.0 

zrms 

i = 

3S.0 

nodhmp = 0 




M0DE = -1 


GAMMA= 

0.0 

NU = 

90.0 

NOE = 


4 

NDREV = 0 




FLARE CONTROL TIME INT 

INCREMENT 

DELAY 









ALPHA 


2.0 

7.0 

0. 

0 




Thrust 4 ej- Vet-Vic.-v 

to 119 s 


POWER 


2.0 

3.0 

1. 

0 





S 4-0 t-hrUST 


TRIM CONFIGURATION 









Ho% e^ecf,/d 


GAMMA= 

-6.40 

ALPHA= 6. 

16 DE 

4.54 

NU= 

90.00 





APPPWR= 

20.0 












AX= 

-0.003 

AZ = 

0.009 









TIME 

DIST 

ALT 

VEL 

AX ( G ) 

AZ<G) 

TH 

TC 

LIFT 

DRAG GAMMA 

ALPHA 

Nil 

DE 

0.0 

0.0 

35.0 

68.0 

-0.00 

0.00 

2240. 

6826. 

4541 0.4 

5535.6 -6.90 

6.16 

90.00 

4.54 

0.2 

22.8 

32.2 

68.0 

-0.00 

0.00 

2240. 

6826. 

45410.4 

5535.6 -6.90 

6.16 

90.00 

4.54 

0.4 

45.6 

29.5 

68.0 

-0.00 

-0.01 

2240. 

6826. 

46039.2 

5772.4 -6.87 

6.86 

90.00 

4.15 

0.6 

68.4 

26.7 

67.9 

-0.01 

-0.03 

2240. 

6826. 

46729.4 

6100.7 -6.79 

7.56 

90.00 

3.74 

0.8 

91.1 

24.0 

67.9 

-0.02 

-0.04 

2240. 

6825. 

47312.5 

6331.3 -6.68 

8.26 

90.00 

3.33 

1.0 

113.9 

21.4 

67.8 

-0.02 

-0.05 

22 39. 

6825. 

47875.1 

6558.5 -6.54 

8.96 

90.00 

2.91 

1.2 

136.6 

18.8 

67.7 

-0.03 

-0.06 

2239. 

6825. 

4R410.1 

6781.1 -6.36 

9.66 

90.00 

2.47 

1.4 

159.3 

16.3 

67.5 

-0.04 

-0.08 

2273. 

6927. 

49070.5 

7063.3 -6.14 

10.36 

90.00 

1.79 

1.6 

181.9 

13.9 

67.3 

-0.05 

-0.09 

2 306. 

7028. 

49643.3 

7369.8 -5.88 

11.06 

90.00 

1.09 

1.8 

204.5 

11.7 

67.1 

-0.06 

-0.10 

2339. 

7129. 

50170.3 

7666.9 -5.58 

11.76 

90.00 

0.38 

2.0 

227.0 

9.5 

66.8 

-0.07 

-0. 12 

2373. 

723 1 . 

50646.9 

7953.1 -5.24 

12.46 

90.00 

-0.36 

2.2 

249.4 

7.5 

66.5 

-0.09 

-0.13 

2406. 

7332. 

51272.9 

8471.1 -4.87 

13.16 

90.00 

-1.13 

2.4 

271.7 

5.7 

66. 1 

-0.09 

-0.12 

2439. 

7433. 

50992.4 

8275.8 -4.47 

13.16 

90.00 

-1.54 

2. 6 

293.9 

4.0 

65.7 

-0. 10 

-0.12 

2472. 

7534. 

50830.7 

8274.0 -4.09 

13.16 

90.00 

-1.98 

2.8 

316.0 

2.5 

65.3 

-0.10 

-0.12 

2 505. 

7634. 

50664.5 

8271.6 -3.72 

13.16 

90.00 

-2.43 

3.0 

337.9 

1.2 

64.9 

-0.11 

-0.11 

2538. 

7735. 

50468.4 

8264.9 -3.36 

13.16 

90.00 

-2.91 

TOUCHDOWN 













3.2 ' 

358.9 

0.0 

64.4 

-0.12 

-0.11 

2571. 

7836. 

60250.3 

8255.2 0.00 

13.16 

90.00 

-3.41 

3.2 

359.7 

0.0 

64.3 

-0.08 

-0.00 

2 504. 

6268. 

35276.0 

3403.1 0.00 

0.00 

94.80 

0.00 

3.4 

381.3 

0.0 

64.0 

-0.08 

-0.00 

2501. 

6211. 

35105.2 

3404.9 0.00 

0.00 

94.98 

0.00 

3.6 

402.9 

0.0 

63. 6 

-0.09 

-0.00 

2434. 

4644. 

30113.3 

3555.4 0.00 

0.00 

99.78 

0.00 

3.8 

424.3 

0.0 

63.3 

-0.10 

-0.00 

2367. 

3077. 

24197. 1 
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5.2 565.3 0.0 

5.6 583.6 0.0 

5.6 601.6 0.0 

5.8 618.7 0.0 

6.0 635.6 0.0 
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7.0 712.9 0.0 
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